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Abstract 

Conditions for a maximum or minimum of Poisson's ratio of anisotropic elastic 
materials are derived. For a uniaxial stress in the 1— direction and Poisson's ra- 
tio V defined by the contraction in the 2— direction, the following three quantities 
vanish at a stationary value: su, \lvs\^ + S25] and [(2i/ — l)si6 + •526]) where s/j 
are the components of the compliance tensor. Analogous conditions for stationary 
values of Young's modulus and the shear modulus are obtained, along with second 
derivatives of the three engineering moduli at the stationary values. The stationary 
conditions and the hessian matrices are presented in forms that are independent of 
the coordinates, which lead to simple search algorithms for extreme values. In each 
case the global extremes can be found by a simple search over the stretch direction 
n only. Simplifications for stretch directions in a plane of orthotropic symmetry are 
also presented, along with numerical examples for the extreme values of the three 
engineering constants in crystals of monoclinic symmetry. 

Keywords: Poisson's ratio, shear modulus. Young's modulus, anisotropy, elasticity 



1 Introduction 

Poisson's ratio z/. Young's modulus E and the shear modulus G, collectively called the 
engineering moduli, are of fixed value in isotropic materials and related by IGiX^v) = E. 
No such connection holds in anisotropic elastic solids, and all they become dependent upon 
the directions of stretch, lateral strain, and the shear directions. Hayes P derived some 
universal relations between values for certain pairs of orthogonal directions. However, 
apart from cubic symmetry j2] there is no general formula for the directions and values 
associated with the largest and smallest values of the engineering moduli. The purpose of 
this paper is to provide systematic methods which can be used to find the extreme values 
of the engineering moduli in any type of anisotropy. 
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The problem of finding the extreme values of Young's modulus is the simplest since 
E depends only on a single direction of stretch. Numerical searching is practical and 
straightforward; thus Cazzani and Rovati provide a detailed analysis of the extrema of 
Young's modulus for cubic and transversely isotropic materials [3] and for materials with 
tetragonal symmetry with extensive illustrative examples. Boulanger and Hayes j3] 
obtained analytic expressions related to extrema of Young's modulus. For stretch in the 
1-direction, they showed that E = 1/su achieves a stationary value if the two conditions 
Si5 = and sie = are satisfied. In a pair of complementary papers, Ting derived 
explicit expressions for the stress directions and the stationary values of Young's modulus 
for triclinic and monoclinic [Hj, orthotropic, tetragonal, trigonal, hexagonal and cubic 
materials [Zj. We will rederive the stationary conditions for E below, along with conditions 
required for a local maximum or minimum. 

Poisson's ratio and the shear modulus depend upon pairs of orthogonal directions, 
which makes their classification far more complicated than for E. At the same time, there 
is considerable interest in anisotropic materials which exhibit negative values of Poisson's 
ratio, also called auxetic materials jHj. In sharp contrast to isotropic solids for which 
— 1 < u < 1/2, the value of u is unrestricted in anisotropic materials and may achieve 
arbitrarily large positive and negative values in the same material. The first hint of this 
surprising possibility was given by Boulanger and Hayes P who presented a theoretical set 
of elastic moduli for a material with orthorhombic symmetry which satisfy the positivity 
requirements, but exhibit simultaneous arbitrarily large positive and negative values of 
Ting and Chen ^U] and Ting subsequently demonstrated that the same remarkable 
phenomenon can be obtained in any non-isotropic material symmetry, including cubic 
symmetry and transverse isotropy. Further explanation of the effect in cubic symmetry is 
provided in Section IHlbelow and in We note that Rovati presented extensive numerical 
examples of auxetic behavior in orthorhombic ^21 and monoclinic materials ^2]; while 
Ting and Barnett ^3] derived general conditions required for the occurrence of negative 
values of u. 

The purpose of this paper is to provide a general framework for finding the maximum 
and minimum values of u and G in anisotropic materials. Some progress in this regard 
is due to Ting [7] who discusses the conditions for extreme values of the shear modulus 
with particular attention to shear in planes of material symmetry. As far as I know, there 
are no results reported to date on conditions necessary for extreme values of the Poisson's 
ratio. Particular attention is given in this paper to the Poisson's ratio, with the emphasis 
on deriving conditions that are independent of the coordinate system used. It will become 
evident that there is a strong analogy between the problems for the shear modulus and the 
Poisson's ratio. In particular, by formulating the problems in a coordinate free manner, 
the task of searching for extreme values of both is similar. 

The outline of the paper is as follows. The three engineering moduli are introduced in 
SectionEl along with the equations for the transformation of elastic moduli under rotation. 
These are used in Section 0] to derive the conditions required for stationary values of u, 
E and G, and the values of the second derivatives (Hessian) at the stationary points 
are determined in Section El The various results are all cast in terms of stretch, strain 
and shearing along coordinate axes. The more general format for stationary conditions in 
general directions, independent of the coordinates, are presented in Sectional The specific 
application to Poisson's ratio is considered in Section IHl The stationary conditions for 
stretch in a plane of orthotropic symmetry are derived, and it is shown that at most 
four stationary values of u can occur, two for in-plane lateral strain, and two out-of-plane. 
These results are applied to the specific case of extreme values of u in materials with cubic 
symmetry, recovering results of [2]. Applications to generally anisotropic materials are 
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also discussed, and a fast procedure for searching for extreme values of v is derived and 
demonstrated for some materials of monoclinic symmetry. Finally, in Section [7| we present 
a similar procedure for finding the global extreme values of G in generally anisotropic 
media, with numerical examples. 

2 Definition of the engineering moduli and prelimi- 
nary equations 

Poisson's ratio measures lateral strain in the presence of uniaxial stress. For any orthonor- 
mal pair of vectors {n, m}, the Poisson's ratio Vnm = i^(n, m) is defined by the ratio of 
the strains in the two directions for a uniaxial state of stress along one of them ^3] : 

£ : mm 

^nm = , forcr = o-nn, (1) 

£ : nn 

where e and cr are the symmetric tensors of strain and stress, respectively, and ab is the 
tensor product, sometimes denoted a ® b. The Young's modulus En = E{n) relates the 
axial strain and stress. 

En = for cr = ann. (2) 

£ : nn 

The third engineering modulus is the shear modulus Gnm = G{n, m), 

Gnm = 7 T for a = a (nm + mn). (3) 

£ : (nm + mn) 

Tensor components are defined relative to the fixed orthonormal basis {ei, e2, 63}, 

The stress 0"^ and strain related by^ 

Sij = Sijki<Jki- (5) 

Here Sijki denote the components of the fourth order compliance tensor. We use the Voigt 
notation for conciseness; comphance is S = [sjj] , /, J = 1, 2, ■ ■ ■ , 6, with 1 = 1,2, 3, 4, 5, 6 
corresponding to ij = 11, 22, 33, 23, 31, 12, and sjj = sjj. 

The goal is to find conditions for a maximum or minimum of the engineering moduli, 
with emphasis on Poisson's ratio. With no loss in generality assume that the n in the ei 
direction, and m is in the e2 direction. Thus, we consider u = 1/12, E = Ei and G = G12, 
i.e.2 

- = -^, i? = -, G = J-. (6) 

Our objective is then to find conditions for a maximum or minimum of each engineering 
modulus under the assumption that the material is assumed to be free to orient in arbitrary 
directions with oriented moduli while the stress remains of fixed orientation. This is 



^Lower case Latin suffices take on the values 1, 2, and 3, and the summation convention on repeated 
indices is assumed unless noted otherwise. 

^Note that we take sge = S1212 although it is common to subsume the factor of 4 in the definition of 
See in eq. ©s- 
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equivalent to stationarity conditions for Unm, En and G„m for a fixed orientation material 
while {n, m} range over all possible orthonormal pairs. 

We therefore need to consider how u, E and G of © vary under general rotation of 
the material. Define the rotation by angle 9 about an arbitrary direction q, |q| = 1 , as 
Q(Q; ^) ^ ^(3), such that vectors (including the basis vectors) transform as r — r' = Qr. 
Under the change of basis associated with Q(q, 0), second order tensors (including stress 



and strain) transform as cr ^ cr' 



where a',,- 



0",- 



Q 



where Q 



ijrs 



■ (^QirQjs ~l~ QisQ jr) • 



(7) 



In order to simplify the algebra we use the connection between fourth order elas- 
ticity tensors in 3 dimensions and second order symmetric tensor of 6 dimensions |15j . 
Accordingly, the 6x6 matrix S with elements Su is defined as 

S = TST, where T = diag(l, 1, 1, ^2, ^2, v^). (8) 

Rotation of second and fourth order tensors is most simply presented in terms of the 6x6 
rotation matrix Q which is the 6-dimensional version of the fourth order tensor Qijrs, 
introduced by Mehrabadi et al. ^Hl- Fourth order tensors transform as S ^ S' = QSQ"'", 
where Q(q, 6') is an orthogonal second order tensor of six dimensions, satisfying QQ^ = 
Q^Q = I =diag(l,l,l,l,l,l). It satisfies 



dQ 

89 



(q,^) = R(q)Q, Q(q,0)=I, 



(9) 



where R is a skew symmetric six dimensional tensor linear in q. 



R(q) 



( ° 











V2q2 


-V2g3\ 











-V2qi 





V2qs 











V2qi 


-V2q2 










-V2qi 





q3 


-q2 







V2q2 


-q3 





qi 


\ V2q3 







q2 


-qi 


/ 



(10) 



Further details can be found in [121 EI . 



3 Extremal conditions 

Consider any one of the engineering moduli, say /, as a function of both the underlying 
compliance and of the rotation Q. A stationary value is obtained if / is unchanged with 
respect to additional small rotations. In order to formulate this more precisely, assume / 
is at a stationary point, and define 



S(q,^) = Q(q,^)SQ^(q,^) 



(11) 
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Define the rotational derivative, 



/'(q) 



de 

ds 



(S(q,^)) 



61=0 



IJ 



(12) 



The elements Sjj(q) of the the rotational derivative of the compliance can be expressed 
by using the representation with (fTTj) . 



S' = R(q)S + SR^(q) 



(13) 



Thus, 







/ 


4si5 


-4sie \ 








-4S24 





4S26 




■S33 




4S34 


-4S35 







^12 




— 2Si4 


2^25 


2Si6 ~ 2^26 








2S24 — 2^34 


— 2^25 


2^536 




■^13 




2Si4 


2^535 — 2Si5 


~2S36 








S12 - Sl3 


-S16 + 2S45 


Si5 — 2S46 




^25 




S26 — 2S45 


■S23 - ■S12 


-S24 + 2S56 








~S35 + 2S46 


S34 — 2S56 


■5l3 ~ ^^23 






s[. 






— Sii + Si3 + 2^55 


— Si4 — 2^56 






"^16 




— Sl5 


Si4 + 2^56 


Sll — S12 — 2Sqq 




(I2 


■5 24 




S22 - S23 - 2S44 


-■S26 


S25 + 2S46 






■^26 




-S25 - 2S46 


■§24 


-■§22 + ■S12 + 2S66 






■^34 




"■533 + -523 + 2^44 


"■■536 2S45 


■535 




■^35 




■536 + 2^45 


■S33 ■5l3 ~ 2^55 


"■534 




S44 




2^24 — 2^34 


— 2^46 


2S45 




"^55 




2S56 


2S35 - 2si5 


-2S45 




■^66 




-2S56 


2S46 


2si6 - 2S26 




■S45 




■S25 - -§35 + S46 


-Si4 + S34 - S56 


S55 — S44 




■^46 




""■§36 + ■526 ~ ^545 


S44 — See 


■5l4 ~ S24 + ■556 




vJ 


\ See — S55 


■^36 ~ ■S16 + ■S45 


"■§25 + ■Sis ~ ^^46 / 





(14) 



Here sjj are the values at ^ = 0, which are independent of q. The derivatives in eq. (jl4j) 
are linear functions of the coordinates of q, and so we may write 

/'(q) = d(^) ■ q, (15) 

where the vector d^'^-' is independent of q and depends only on the compliances. 

The engineering modulus / is stationary with respect to the direction n, and the 
direction m where applicable, if /'(q) vanishes for all axes of rotation q. This is equivalent 
to requiring that all possible deviations in n and m leave / unchanged to first order in 
the rotation. The stipulation that this hold for all rotation axes covers all permissible 
transformations. The general condition for stationarity is therefore that the vector d*^-^) 
must vanish, i.e., 

d^-'^-' =0 at a stationary point of /. (16) 

We now apply this formalism to the three engineering moduli and derive d'^''^ d''^^ and 
d*^'^'' in turn. 
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3.1 Poisson's ratio 

For Poisson's ratio, eq. (jl2p becomes 

'^'(q) = [si2 s'niq) - Sn s'^^i^l)]- (17) 
Thus, from eqs. (fT^ through (fTTj) . we have 

d^'') = 2s]f [siiSuei + (2si2Si5 - SiiS25)e2 + (siiS26 - S11S16 - 2si2Si6)e3] . (18) 

Setting this to zero and using the definition of z/ in and the fact that sn > 0, we 
obtain three conditions for a stationary value of Poisson's ratio: 

Si4 = 0, (19a) 
2z/si5 + S25 = 0, (19b) 
(2z/ - l)si6 + S26 = 0. (19c) 

These must be simultaneously satisfied at a maximum or minimum of u. Note that the 
compliance elements appearing in ()19|) are all identically zero in isotropic materials. 

3.2 Young's modulus 

Proceeding in the same manner as for the Poisson's ratio, and using E'{q) = —Siis'ii{q), 
gives 

d(^) = 4s^^ ( - si5e2 + sieea) . (20) 
Setting this to zero implies the conditions for an extremum in Young's modulus 

si5 = 0, (21a) 
sw = 0. (21b) 

These agree with two conditions determined by Boulanger and Hayes [S] and by Ting 

3.3 Shear modulus 

The rotational derivative of the shear modulus is G"(q) = — |s^g^Sgg(q), and the gradient 
vector is 

d^"^^ = -SqI [Ssgei - 54662 + (^26 - 516)63] . (22) 

Hence, the shear modulus has an extreme value if the following three conditions hold: 

S56 = 0, (23a) 
Sag = 0, (23b) 
SI6 - 526 = 0. (23c) 
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4 Second derivatives 

The nature of a stationary value of the general engineering modulus / can be discerned, at 
least locally, by the second derivative. By analogy with eq. fll2|) . we define the rotational 
second derivative. 

The elements s'/j(q) of the the rotational second derivative of the compliance follow from 

S"(q) = R^S + SR2^ + 2RSR^ . (25) 

This is a direct consequence of eq. (jl3|) . We do not need all 21 elements, and for 
brevity only present the following three values which are necessary to evaluate the second 
derivatives of the engineering moduli, 

s'/i = 4[(si3 - -511 + 2s55)g^ + (si2 - Sii + 2s66)g3 ^ 2(si4 + 2s^&)q2qz + sisgs^i + ■5i6gig2] , 

(26a) 

s'[2 = 2[(Si3 - Si2)ql + (S23 - 512)^2 + (^ll + S22 - 2Si2 - 4s66)g| 

+ (Si4 - 2S24 + 4s56)g2g3 + (^25 - 2Si5 + 4546)^3?! + (^16 + ^26 " 4s45)gig2] , 

(26b) 

See = 2[(s55 - see)?! + (-544 - S6e)g2 + (-^le + S2e - 2s3e - 2545)^1^2 
+ {su + S22 - 2si2 - 4se6)g3 + (2s25 - 2si5 + 3s4e)gig3 

+ (2si4 - 2s24 + 3s5e)g2g3] • (26c) 

Before applying these to the three engineering moduli f = u, E and G, we note that 
in each case that / is a homogeneous function of degree or - in the compliance elements. 
Consequently the second derivative /" evaluated at the stationary point where /' = 
simplifies because the first term in ()24j) vanishes, leaving 



n^) = -^s'M at/'(q) = 0. (27) 

The terms s'/j(q) are second order in q, and we can write 

/"(q) = D(/) : qq at /'(q) = 0, (28) 

where D^-^^ = D*^-^^"^ is a non-dimensional symmetric 3x3 matrix which is independent 
of q. Thus, D*^-^) is positive (negative) semi-definite at a local minimum (maximum) of /. 
The condition for a local minimum (maximum) is therefore that the three eigenvalues of 
D^'^^ are positive (negative). If the matrix is not definite and has eigenvalues of opposite 
sign, then the modulus has a locally saddle shaped behavior. 

4.1 Poisson's ratio 

The second derivative of the Poisson's ratio at a stationary point is, using ()27|). 

z/"(q) = -sr/ {s'l^ + lys';,) at u'{q) = 0. (29) 
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Thus, when i^'(q) = 0, eqs. and give 

'^"(q) = "1- I ■5ii(si2 - Si3)ql + [■5ii(si2 - Sis) + 2Si2(Si3 - Sn + 2s55)]g^ 
^11 I 

+ [sii(2si2 + 4s66 - Sn - S22) + 2si2(si2 - Sii + 2sQe)]qi 
+ 2[sii(s24 - 2S56 - ^Su) - 2Si2(Si4 + 2s.,Q)]q2q3 

+ 2[sii(si5 — 2^46 — 2^25) + SuSi^^q^qi 

+ 2[sii(2s45 - ^si6 - ^sse) + si2Si6]gig2| • 

Since this is evaluated at the stationary value, we may use to simplify and obtain 

/ S12 - Si3 2S45 - S16 Si5 - 2S46 \ 



(30) 



2S45 ~ S16 



S12 

2i^(si3 - 



- S13- 



2S55) S24 - 2(1 - 2Z/)S56 



\si5 - 2s46 S24 - 2(1 - 2z/)s 



S12 - S22 + 2S66 + 
56 (1 - 2z/) (Si2 - Sn + 2S66) / 



(31) 



where we have used the definition of u to simplify the elements. 

4.2 Young's modulus 

At a stationary point 

E"{q) = -s^^ s'/i at E'{q) = 0. 
Using eq. (j26p and the extremal conditions (|2ip. we find 

^0 



'11 



Sn - Si3 - 2S55 Si4 + 2S56 

,0 Si4 + 2S56 Sn-Si2-2S66, 



(32) 



(33) 



Note that D*^-^) is rank deficient (of rank 2), which is a consequence of the fact that E 
is invariant under rotation about the ei stretch axis. The local nature of the stationary 
value depends upon the two non-zero eigenvalues of the matrix. 



4.3 Shear modulus 

The shear modulus satisfies 



and hence. 



2^2 



G"iq) = -^s^e'4, at G'(q) = 0, 



See ~ S55 S36 — S16 + S45 Si5 — S25 

Sse ~ Sie + S45 See ~ S44 S24 — S14 



(34) 



(35) 



Sl5 - S25 



S24 - Si4 2si2 + 4s6e - Sn - S22, 
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5 Coordinate invariant formulation 

In this section we rephrase the results for the stationary conditions and for the second 
derivatives at the stationary conditions in coordinate invariant form. Let {n, m, p} be an 
orthonormal triad analogous to {ei, 62, 63} before. Define the non-dimensional symmetric 
second order tensors A, B, N, M and P, as follows: 

^ij ~ ^nn^ijkl^k^h B^j = Sj^^SijklTni^mi, Cij = S^^SijklPkPh (36a) 

^ij = s~^Sikjinkni, Mij = s'^Sikjiirikmi, Pij = s'^SikjiPkPi- (36b) 
where s„„, = Si^kiniUjUkni. Thus, 

i^(n, m) = — A : mm, -E'(n) = 1/ s„„, ^(n, m) = -E(n) / (4N : mm) . (37) 

5.1 Poisson's ratio 

The derivative of Poisson's ratio can be expressed in general form as 

d*^'') = [A : (mp+pm)]n-[(2z/A+B) : (pn+np)] m+{ [(2z/-1)A+B] : (nm+mn)}p, 

(38) 

which follows from eq. (fT^ . and may be checked by the substitution {n, m, p} 
{61,62,63}. This provides a local expansion of the Poisson's ratio for small rotation 
about the axis q: 

z/(q, 9) = z/(q, 0) + d^'^) ■ q ^ + 0(9^) . (39) 

In particular, d'^'^^ = at a stationary point. 

The second order tensor D^^^ of ()31|) becomes, in general format, 

D*^'") = [A : (mm - pp)]nn + [u - [(1 + 2i^)A + 4z/N] : ppjmm 
+ {2u^ - 1 + [4(1 - z/)N - M] : mm}pp 

+ [(2P — N) : (nm + mn)] ^(nm + mn) 
+ [(N - 2M) : (pn + np)] ^(pn + np) 

+ { [M - 2(1 - 2u)N] : (mp + pm)}^(mp + pm) . (40) 
Hence, the local expansion near a stationary point is 

i^{q, 9) = z/(q, 0) + io^'^) : qq ^2 ^ 0{9'') . (41) 

5.2 Young's modulus 

In the same way as before we find that 

^(E) ^ 2E{n) [ - A : (pn + np) m + B : (nm + mn) p] , (42) 

and 

D^^^ = -E'(n) { [a : (nn — pp) — 2P : nn] mm + [A : (nn — mm) — 2M : nn] pp 

+ (A + 2N) : (mp + pm)}^(mp + pm)}. (43) 
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5.3 Shear modulus 



Similarly, for the shear modulus 



d(^) = 4 



[N : (mp + pm) n — M : (pn + np) m + (B — A) : (nm + mn) p] , (44) 



E{n) 



and 



8 



B — — — {N : (mm — pp) nn + M : (nn — pp) mm + 4N : mm pp 

E{n) ^ 

+ (P + C — A) : (nm + mn)-(nm + mn) + (A — B) : [(mm — nn) pp 

1 1 
- (mp + pm)-(mp + pm) + -(pn + np)-(pn + np)] }. 




(45) 



6 Applications to Poisson's ratio 

We now concentrate on general properties of the Poisson's ratio, applying the formalism 
for the stationary value to different situations. We begin by consider the general case of 
a plane of material symmetry in an orthotropic material. 



6.1 Plane of symmetry in orthotropic material 

We assume the stretch direction n lies in a plane of symmetry of an orthotropic material, 
and the direction of contraction m lies (a) perpendicular to the plane or (b) in the plane. 
This configuration includes all planes in hexagonal materials that contain the axis of 
symmetry, and therefore provides the stationary values of u in materials with hexagonal 
symmetry (transverse isotropy). 

With no loss in generality, let e^^^ be the normal to the plane of symmetry and let 
case (a) correspond to m = 63 and case (b) corresponds to m in the plane of e^°\ In 
both (a) and (b) n lies in the plane of e^i\e'^\ Define the rotated axes 



ei = cos 6* e^^ + sin 6* 62 , e2 = — sm9 e{ + cosO 62 , 63 = eg . (46) 



Let Sjj denote the compliances relative to the fixed set of axes {e[ ,63 }, and s/j 
the compliances in the coordinates of the rotated axes. By definition of a symmetry plane. 
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all coefficients Sijki with index 3 appearing once or thrice are zero. Then, 

^11 = 0,0 + 7 (-^ii + •^22 - 25o) [ ^ - COS 2e\ , (47a) 

Oil + 022 ~ -^OO ^ oil + 022 "~ 

512 = S12 + ^{Su + ^22 - 2So) sin2 2^, (47b) 

513 = ^ (5l3 + ^23) - ^ (^23 - ^13) COS 2^, (47c) 

sie = ]iSn + S22-2So) L ^'V^"oc " 2^' (47d) 

4 Du + D22 — ^00 

526 = 7 (^11 + ^22 - 25o) L ^'V^"oc + 20, (47e) 

4 Oil + 022 — ^00 

536 = ^(^23 -^13) Sin 20, (47f) 

where 

So = S,2 + 2See. (48) 
Thus, in the two cases to be considered, we have n = rii, and sn = 1/E{6). 

(a) m = 63 perpendicular to the plane of symmetry 

We now consider stationary values of u = u^^, for which the three conditions for stationary 
i> are, instead of (|TTH) . 

Si4 = 0, 2l^Si6 + S36 = 0, {2U - l)si5 + S35 = . (49) 

The ffist and last of these are automatically satisfied, based on the assumed material 
symmetry. Using Sig and s^q from (j37j), the second of ()49p implies that sin 20 = 0, which 
is the exceptional case of prior symmetry, or that 6 satisfies 

cos 20 = ~ + -{ ^^3-^13 ^ 

5*11 + 5*22 — 25*0 ly Sn + 6*22 — 25*0 

Using this to eliminate cos 20 from the expression for u = = —Srs/su yields a quadratic 
equation for z/, 

i/^ - - -pa = 0, (51) 

where 

_ ('S'o — 'S'22)>S'i3 + (5*0 — 5'ii)5'23 _ (5*23 — 6*13) 

^a. — Q 3^ 7;^ ; Pa — ~^ Q 7^2" ■ \^^) 

OIIO22 — Oq O11O22 — Oq 

Define E* and 0* by 

'S'll + ^22 — 25*0 



E* = "'7^ ^2 , (53) 

011O22 — Oq 

then we may identify E* as the value of E{6) when the second term on the RHS of ()47a|) 
vanishes, i.e. E* = E{6*), where 0* satisfies 
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The angle 9* defines the direction at which E is stationary (maximum or minimum). It 
exists iff the RHS of ()54|) hes between —1 and 1. Regardless of whether or not the angle 
exists, it can be checked, Ua = i^rsi^^*) = —E*Si3{9*). The value of Young's modulus in 
the stretch direction 6 for the stationary value of u satisfies 

^ + ^=2. (55) 

In summary, the possible stationary values for stretch in the plane of symmetry and 
the strain measured in the direction perpendicular to the plane are 

Va± = \va±\[vl^Pafl\ (56) 

and the stationary values occur if — 1 < ^a± < 1, where 

^ _ S22 — Sii ^ 1 ^ '5'23 ~ '5'l3 ^g^^ 

5*11 + 5*22 — 25*0 I^a± 'S'li + 5*22 — 25*0 

in which case the direction of stretch is given by 6' = | cos~^ 7a±. Otherwise the stationary 
values occur at 6* = and n/2. 

(b) m = 62 in the plane of symmetry 

With n = ei again, conditions ()19a|) and (|19b|l are met and only ()19c|) is not automatically 
satisfied. Substituting the expressions for sie and S26 from (jTfjl into equation ()19c|l implies 
that either sin 26 = 0, which is simply the axial case, or 

^22 — i^ii 1 / 5*22 — "S'li 



COS 29 



+ [ ~ ) (5^ 

u — 1 Sii + 5*22 ^ 25*0 



•S'li + 5*22 — 25*0 

Substituting this into the relation for u = z/12: lysnn + S12 = and eliminating cos 20 
produces a quadratic equation in u. The equation is most simply expressed as a quadratic 
in the shifted Poisson's ratio (i^ — 1): 

{u-iy-{uh-l){u-l)-^Ph = 0, (59) 

where 

H = -l-E*{Sr,-So), P6= IfV^il - (60) 

'-'ll'-'22 — ^0 

That is, z/fe = uuiO*) = —E*si2{0*). Note that the value of Young's modulus E in the 
stretch direction 6 satisfies 

E* Ub-l 

+ J^ = 2. 61 
E u — 1 

In summary, the possible stationary values for stretch and strain both in the plane of 
symmetry are 

i^b± = ^iH + l)± ^{{u, - 1)2 + p,y/\ (62) 
and the stationary values occur if —1 < 7^,-1- < 1, where 



7b± 



5*22 ~ "S'li 1 / 5*22 — S 



Sii + 5*22 — 25*0 i'b± — 1 5*11 + 5*22 — 25*0 



^ 022 - on ^ ^gg^ 
1 S'li + S')'?. ~ 2Sr) 



in which case the direction of stretch is given hj 9 = ^ cos ^ 7b±. 
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6.2 Example: cubic materials 

The general coordinate invariant form of the comphance of a cubic material is 

S = (— - — ) J + — I + ^) D, 



(64) 



where I and J are fourth order isotropic tensors, lijki = |(5ifc5jz + ^u^jk)-, Jijki = k^ij^ki, 



and 



D = aaaa + bbbb + cccc . 



(65) 



Here, the orthonormal triad {a, b, c} is coaxial with the cube axes. The condition that 
the elastic strain energy is always positive definite is that the three moduli k, Hi and 
yU2 are positive. We use, for simplicity, crystallographic-type notation for unit vectors, 
e.g., pqf where p, q and r are positive numbers, indicates the unit vector (pa + gb — 
rc) / + + r'^ . 




Figure 1: The irreducible l/48th of the cube surface is defined by the isosceles triangle with 
vertices as shown. 



We consider the stationary conditions ()19|) for stretch in the ei direction and lateral 
strain in the 62 direction. Since S14, S15, sie, S25 and S26 vanish for isotropic media, it 
follows that the only contribution to these quantities is from the tensor D. We may 
therefore rewrite the stationary conditions ()19p as 

Di4 = 0, (66a) 
2z/Di5 + D25 = 0, (66b) 
{2V - l)Di6 - ^26 = 0, (66c) 

where Dn = -D1123, etc. The realm of stretch directions that needs to be considered may 
be reduced to those defining the irreducible l/48th of the surface of the unit sphere. This 
in turn is defined by l/48th of the surface of the cube, see Fig. ^ where the vertices of the 
triangle correspond to the directions 001, 110, and 111. In a separate paper |2I it is shown 
that the extreme values of u do not occur within the interior of the triangle. It turns out 
that the extreme values are only possible for stretch direction ei along 001, 110, or in 
certain cases, for ei located along the edge between 110 and 111. In the case that ei = 001, 
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the lateral direction 62 may be any orthogonal direction, and when ei = 110 the lateral 
directions are 001 or 110, each of which can correspond to the minimum or maximum for 
u, depending on the elastic parameters k, /ii and fi2- It is clear from the symmetry of the 
situation that the quantities D14, Z^is, Diq, D25 and D2Q vanish identically for ei along 
001 or 110 with 62 as described. A full description of the possible extreme values of u in 
cubic materials is involved but complete, and we refer to |2] for details. 

To summarize the findings of Norris |2] regarding solutions of eqs. ()66|) : all possible 
stretch directions which solve the three stationary conditions are confined to the symmetry 
plane with normal 110, and equivalent planes of symmetry. We can therefore apply 
the results for a plane of orthotropic symmetry. Explicit calculation from ()64|) gives 
^66 = (4/ii)-\ ^12 = ^23 = (9«:)-^ - (6/i2)-\ 5i3 =_Si2 + x/4, So = 5i2 + 2^66, 
-^ii = 5'o + x/4, 5'22 = So + x/2, where x = {f^2) ^ - (l^i) ^- Thus, 

l^a = l^b = l^Ul = ——, Pa = Pb = P= + 1), (67) 

6k + 2/ii 6 ^p2 

where Um = z/(lll,m) is the Poisson's ratio for stretch in the 111 direction, and is 
independent of the lateral direction m. 

The actual values of the possible extrema for u can be obtained from equations and 
fl^Uj) . (jU^ and (jnZj). Skipping over the unedifying details, see it can be shown that 
only the stationary values and z/b_|_ are possible global extrema, 

1 1 



-^^111 - -y^^iii +P, (68a) 
I (^111 + 1) + ^ ^J{uni-lY + P, (68b) 



Note that um is independent of /i2, which only enters these expressions via the term p. 
The extremely large values of Poisson's ratio discovered by Ting and Chen ^U] correspond 
to p 3> 1, which can occur if P2I Pi ^ 1- Under this circumstance Ua- is large and negative, 
1^5+ is large and positive, and the magnitudes are, in principle, unbounded [TUl El 12] • 
The directions associated with the global extrema are given by 

cos20,_ = ^--^, cos2e,+ = \+ ^ (69) 
3 3z/„_ 3 3(z/6+ - 1) 

Both directions bifurcate from 6' = [2j, and therefore these extreme values only occur 
ifi^<— |orz/>|, respectively. A complete description of the extrema for all possible 
values of the elastic moduli is given by Norris [2j. 



6.3 Application to generally anisotropic materials 

We first present a result which suggests a simple algorithm for searching for global extreme 
values of v in generally anisotropic materials. 

6.3.1 A local min-max result for Poisson's ratio 

The tensor of second derivatives of z/(n, m) at the stationary point, D*^'^^ must be posi- 
tive or negative definite in order that the stationary point is a minimum or a maximum. 
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respectively. Consider a possible minimum, then a necessary although not sufficient con- 
dition is that the three diagonal elements of D'^'^) are positive. In particular, eq. (pn|) 
gives 

D^'^) : nn = v^p - Vnm > 0. (70) 

This implies that u = Unm must be strictly less than i/„p. At the same time, the stationary 
condition d^*^^ = must hold, and in particular, 

d*^'^^ • n = =^ £ : mp = for cr = a nn. (71) 

This implies that the shear strain Smp in the m, p-plane is zero. Hence, for any unit vector 
r ± n, we have 

£ : rr = (r ■ m)^ e„^m + (r ■ p)^ (72) 

or 

Unr = (r ■ m)^ Unm + (r ■ p)^ Unp 
= I'nm + (r • p)^ [Unp " I'nm) 

> ^nm, (73) 

with equality only for r = m. Thus, we have: 



Lemma 1 Ifvnm is a minimum (maximum) value, then it is also a minimum (maximum) 
among all possible Vnr for r in the plane perpendicular to n. 

This result is a direct consequence of the general expression for 0*^^^ It implies that 
if we can satisfy (fTTj) then the values of and Unp are the extreme values for the given 
stretch direction n. We next show how this single condition can be achieved. 

6.3.2 Satisfaction of one extremal condition 

The stationary values of Poisson's ratio occur, in general, for stretch directions at which 
the vector d*^^^ of eq. ()H8j) vanishes. We now show that one of the three components can 
be made to vanish; specifically, d'^'^^ ■ n = A : (mp + pm) is zero for an appropriate choice 
of the orthogonal directions m and p. 

We use the fact that the appropriate pair m and p correspond to stationary values of 
A : mm and A : pp. These satisfy A : mm + A : pp =trA — 1 so a maximum in one 
implies a minimum for the other. In order to find these directions for a given n, consider 
the function of m: 

g{m) = A : mm — Am ■ m — 2am ■ n . (74) 

Setting to zero the gradient with respect to m implies that m satisfies 

m = a(A-Al)^^n. (75) 

The scalars A and a follow by requiring that m ■ n = and m ■ m = 1, respectively. The 
former implies that A satisfies 



n ■ (A - AI) n = . (76) 
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This condition can be rewritten by expanding the inverse in terms of the cofactor matrix 
of (A — Al), and using the property A : nn = 1, which yields a quadratic in A: 

A2 + A(l -trA) +adj(A) : nn = 0. (77) 

Here adj(A) is the adjoint matrix. If A is invertible, then adj(A) = (detA)A~^, and 
more generally the adjoint is the transpose of the cofactor matrix. Thus for any type of 
anisotropy, and for any direction n, it is a straightforward to determine the appropriate 
orthogonal pair m, p, which automatically give d^*^^ ■ n = 0. Finding the right axes 
requires solving the quadratic (fTTjl for A. The associated values of Poisson's ratio are the 
maximum and minimum for the given direction n. In this way, finding extremal values of 
Poisson's ratio for all possible n is reduced to seeking values which satisfy the remaining 
two conditions: 

(2z/A + B) : np = 0, [{2iy - 1)A + B] : nm = 0, (78) 

or equivalently, 2z/si5 + S25 = and (2i/ — I)si6 + S26 = 0, respectively. This is the strategy 
used to determine the global extrema of u in materials with cubic symmetry j2]. 

6.3.3 Algorithm for finding global extreme values of v 

Rather than searching for directions n which satisfy the three stationary conditions on u, 
Lemma ^ suggests that a simple search for maximum and minimum values of Poisson's 
ratio can be effected as follows. For a given n define the pair i^±(n) by 

z/^(n) = -A : m''^^m''^\ (79a) 
mi^) = II (A - Ai^)l)"'n|ri (A - Ai^^l)"'n, (79b) 

Ai^^ = ^(trA-1) ±^[(trA-l)'-4adj(A) :nn]'/'. (79c) 

The search for global extrema is then a matter of finding the largest and smallest values of 
iy±{ii) by searching over all possible directions n. In practice, even for triclinic materials 
with no symmetry, the search only has to be performed over half of the unit sphere 
||n|| = 1, such as n ■ e > for some fixed direction e. In this way, the numerical search is 
equivalent in complexity to that of finding the global extrema of Young's modulus. 

This algorithm was applied to data for two crystals of monoclinic symmetry: Cesium 
dihydrogen phosphate and Lanthanum niobate, with the results in Table 1. The numerical 
results were obtained by using a 100x100 mesh for the hemisphere ||n|| = 1, n ■ e2 > 0. 
It was found that the extreme values of the engineering moduli are not sensitive to the 
mesh size, although the values of the directions do change slightly. 

7 Applications to the shear modulus 

If we compare the vector derivative function for the shear modulus, d^*^) of with 
d*^*^) of (jSHI), we note that the components in the n direction have a similar form, but 
with different matrices involved. Thus, it is N for d*^"-^) ■ n, while for d*^"^) ■ n it is A. 
The properties of the second derivatives D^*^) : nn and D'^'"^ : nn are similarly related. 
Proceeding with the same arguments as for the Poisson's ratio, we may deduce 
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(LaNb04) 


E • 

^mm 


2.27 


-0.50 


0.00 


0.87 




^max 


72.57 


-0.28 


0.71 


-0.64 




^min 


0.70 


-0.96 


0.00 


0.27 




^max 


3.95 


0.00 


1.00 


0.00 




^min 


-3.01 


0.00 


1.00 


0.00 



Table 1: Extreme values of E^ G and v for two materials of monoclinic symmetry with 
symmetry plane n2 = 0. Units of E and G are GPa. Cesium dihydrogen phosphate: Sn 
= 1820.0, S22 = 103.0, S33 = 772.0, S44 = 33.25, S55 = 112.5, sge = 29.25, S12 = -219.0, 
si3 = -1170.0, S23 = 138.0, si5 = 124.5, S25 = -75.0, S35 = -90.5, s^q = 8.25. Lanthanum 
niobate: Su = 66.8, S22 = 14.8, S33 = 146.0, S44 = 5.7, S55 = 265.0, see = 4.675, su = 
16.9, si3 = -94.8, S23 = -30.8, S15 = 118.0, S25 = 45.6, S35 = -186.5, S46 = 0.95. Units in 
(TPa)~^ (Data from Landolt & Bornstein [TU] ) 



Lemma 2 IfGnm is a minimum (maximum) value, then it is also a minimum (maximum) 
among all possible Gnr for r in the plane perpendicular to n. 

This in turn leads to a similar method for finding the global extrema of G. 

7.1 Algorithm for finding global extreme values of G 

Define the pair ^^^(n) by 

^±(11) = ^(n)/(N : mf ^), (80a) 
) = II (N - )l)-V|ri (N - )l)-'n, (80b) 

A^"^) = -(trN-1) ±-[(trN-l)'-4adj(N) :nn]'/'. (80c) 

Then the global extreme values of G can be found by searching over n only. Thus, 
the problem of finding the global extrema of the shear modulus is reduced to the same 
level of complexity as searching for the maximum and minimum Young's modulus. Table 
1 summarizes the extreme values of G found using this algorithm for two crystals of 
monoclinic symmetry. 
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8 Conclusions 

The results of this paper provide a consistent framework for determining the extreme 
values of the three engineering moduli in materials of any crystal symmetry or none. 
General conditions have been derived which must be satisfied at stationary values of z/, E 
and G. These are equations ()19|) . (j2U) and ()23|) . which have also been cast in forms that 
are independent of the coordinates used, in equations (jHHj) . (jl^ and (HH), respectively. 
The associated three hessian matrices which determine the local nature of the stationary 
value, maximum, minimum or saddle, are given in equations ()40|) . PHj) and (jl^ . The 
stationary conditions for Poisson's ratio simplify for stretch in a plane of orthotropic 
symmetry, for which there are at most 4 stationary values of u. Two can occur for in- 
plane stretch and strain, and the other two for out-of-plane strain. This implies that 
transversely isotropic materials have at most four stationary values of Poisson's ratio. 
The results for the plane of symmetry also reproduce known results for cubic materials 
The hessian matrices for u and G lead to algorithms for finding the extreme values. 
The key is to remove the dependence on the m direction by explicit representation of the 
maximum and minimum for a given n direction. The algorithms have been demonstrated 
by application to materials of low symmetry. 
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